Nonlinear control of chaotic walking of atoms in an optical lattice 
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Abstract. - Centre-of-mass atomic motion in an optical lattice is shown to be near the resonance a 
chaotic walking due to the interplay between coherent internal atomic dynamics and spontaneous 
emission. Statistical properties of chaotic atomic motion can be controlled by the single parameter, 
the detuning between the atomic transition frequency and the laser frequency. We derive a Fokker- 
Planck equation in the energetic space to describe the atomic transport near the resonance and 
demonstrate numerically how to manipulate the atomic motion varying the detuning. 



Introduction. — Let us consider a cold two-level 
atom in a laser standing wave forming a ID optical lat- 
tice. Depending on the values of atomic and lattice pa- 
rameters, different regimes of the centre-of-mass motion 
can be identified: oscillations in wells of the optical po- 
tential, ballistic flights with acceleration and deceleration, 
velocity grouping, Brownian motion and random walking 
with Levy flights [1-10]. Theoretical study and experi- 
mental realization of these effects have been done, mainly, 
in the context of laser cooling of atoms when the laser 
is far detuned from the atomic transition frequency and 
the excited state can be adiabatically eliminated. Near 
the resonance, coherent interaction of atoms with a laser 
field may strongly affect the atomic motion. It has been 
shown in Rcfs. [11-13] that, when neglecting spontaneous 
emission (SE), there should exist a deterministic chaotic 
transport of atoms with exponential sensitivity to small 
variations in initial conditions and/or lattice parameters. 
In our recent work [14] we developed a semiclassical the- 
ory of this phenomenon resembling a random walking but 
without any random forces and noise. 

In this Letter we use a Monte Carlo wavefunction ap- 
proach to take into account SE and demonstrate a new 
type of motion in an optical lattice, a chaotic atomic 
walking, with the properties both of the chaotic atomic 
transport, caused by coherent atomic dynamics, and of 
a random process due to SE kicks. It is a kind of ran- 
dom walking with specific statistical properties that can- 
not be classified neither as deterministic chaotic motion 



nor as normal diffusion nor as sub(super)diffusion and nor 
as Levy flights. We derive a Fokker-Planck equation in 
the energetic space, study the statistical properties of the 
chaotic walking and demonstrate numerically how to ma- 
nipulate these properties varying the atom-field detuning. 
Small changes in this parameter may affect drastically the 
atomic transport, transforming the atomic motion from a 
practically regular one to anomalous diffusion. 

Monte Carlo wavefunction approach. We start 
with the non-Hcrmitian Hamiltonian of a two-level atom 
interacting with a strong standing-wave ID laser field 
forming an optical lattice. In the frame rotating with the 
laser frequency w/ , it has the form 
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where a±^ z are the Pauli operators for the internal atomic 
degrees of freedom, X and P are the atomic position and 
momentum operators, uj a , ujf and f2 are the atomic tran- 
sition, laser and Rabi frequencies, respectively, and V is 
the spontaneous decay rate. Internal atomic states are de- 
scribed by the wavefunction = a(t)\2) + b(t)\l), with 
a and b being the complex-valued probability amplitudes 
to find an atom in the excited \2) and ground |1) states. 



Note that the norm of the wavefunction, 



+ \br, is not 



conserved due to non-Hcrmitian term in the Hamiltonian. 
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To study the atomic dynamics in the optical lattice 
we use the standard Monte Carlo wavefunction technique 
to get the most probabilistic outcome that can be com- 
pared directly with corresponding experimental observa- 
tions with single atoms. The method is based on the evo- 
lution of an atomic state \^(t)) while a continuous mea- 
surement of radiation-field state is performed by an ideal 
photodetcctor. The evolution consists of two parts: (1) 
jumps to the ground state (a = 0, \b\ 2 = 1) each of which 
is accompanied by the emission of an observable photon at 
random time moments with the mean time (|a| 2 T) -1 and 
(2) coherent evolution with continuously decaying norm of 
the atomic state vector without the emission of an observ- 
able photon. The decay of the norm of the state vector 
is equal to the probability of spontaneous emission of the 
next photon. This coherent decay without emission of a 
photon is usually interpreted within the context of the 
quantum measurement theory. Let us introduce the new 
real-valued variables (normalized all the time) instead of 
the amplitudes a and b (renormalizcd only after SE events) 
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which have the meaning of synphase and quadrature com- 
ponents of the atomic electric dipole moment and the pop- 
ulation inversion, respectively. We stress that the length 
of the Bloch vector, u 2 + v 2 + z 2 = 1, is conserved. 

The average atomic momentum is supposed to be large 
as compared with the photon momentum Kkf, so the 
translational atomic motion can be treated classically us- 
ing the Hamilton equations. The whole atomic dynamics 
is governed by the following Hamilton-Schrodinger equa- 
tions [15]: 

oo 

X = U) r p, p = —U Sin X + Pj 5(t — Tj), 

3=1 

oo 



u = Av H — uz — uy S(t — Tj), 

3=1 
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where x = kf(X) and p = (P)/Kkf are the normalized 
atomic centre-of-mass position and momentum, respec- 
tively. The dot denotes differentiation with respect to 
the normalized time r = fit. The values of the nor- 
malized decay rate 7 = T/D, and the recoil frequency 
uj r = hk 2 /m a il <C 1 are fixed in this paper, and the nor- 
malized detuning A = (ujf — w a )/f2 is a variable param- 
eter. In the equations of motion (3) Tj are random time 
moments of SE events and pj are random recoil momenta 
with the values between ±1 in a one-dimensional case. 



In terms of the normalized time r the mean frequency of 
SE events is equal to 7(2: + l)/2. At the time moments 
t — Tj, the atomic variables change as follows: p — > p+pj, 
u — > 0, v ^ 0, z — > — 1. The well-known effects of ac- 
celeration, deceleration and velocity grouping have been 
successfully simulated with Eqs. (3) in Rcf. [15]. In the 
present numerical simulation we take a cesium atom with 
the working transition 65*1/2 — 6P3/2 (A = 852.1 nm and 
r = 3.2 • 10 7 Hz) interacting with a rather strong field 
with the resonant Rabi frequency SI = 10 10 Hz. Thus, the 
corresponding normalized recoil frequency is uj r = 10 -5 
and the spontaneous decay rate is 7 = 3.3 • 10~ 3 . 

It follows from Eqs. (3) that the centre-of-mass motion 
is described by the equation for a nonlinear physical pen- 
dulum with a frequency modulation caused by coherent 
internal atomic dynamics and random kicks of the mo- 
mentum. Besides SE recoils any atomic trajectory is de- 
termined by the coherent evolution of the synphase com- 
ponent of the electric dipole moment u between the events 
of SE and its jumps at random time moments Tj. These 
jumps are not small. In our recent paper [14] we devel- 
oped a theory of atomic transport in an optical lattice 
(in the absence of SE) based on a specific behavior of the 
variable u which performs shallow and fast oscillations be- 
tween the nodes of the standing laser wave and changes 
suddenly its value when atoms cross the nodes. The theory 
predicts deterministic chaotic transport at small values of 
the detuning |A| <C 1 whose statistical properties are very 
well described by a stochastic map for the deterministic 
variable u [14]. In fact, atom moves in a rigid optical lat- 
tice just like as in a random optical potential with a com- 
plicated alternation of oscillations in potential wells and 
flights over many wells when it can change its direction 
of motion many times. SE causes further complication of 
this motion. 

Chaotic walking at small detunings, A <C 1. 

Near the resonance (|A| <c 1 ), the following quantity is 
almost conserved between any two acts of SE: 
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where Tj ^ r < Tj + \. H is the total atomic energy which 
is a constant in the purely Hamiltonian system, i.e. with- 
out any relaxation [14]. The last term in (4) with the 
averaging over a time exceeding the period of the Rabi 
oscillations compensates the relaxation. The energy H 
changes suddenly at the moments of SE and decays lin- 
early in between (Fig. 1) whereas the pseudocnergy H 
changes suddenly as well at the same moments but is ap- 
proximately a constant in between. At H > \u\, the atom 
moves ballistically, and its momentum cannot be zero be- 
cause at small detunings |A| the kinetic energy is larger 
than all the other terms in (4). If H < 0, the atom changes 
surely the sign of the momentum during its motion. Thus, 
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Fig. 1: Time evolution of the atomic energy H (dashed line) 
and the synphase component of the electric dipole moment uO 
(solid line) at the same small value of the detuning A = —0.001 
but with different initial conditions. 



if we could construct a mapping for the pseudoenergy Hj , 
we would know approximately when atoms move ballisti- 
cally and when they turn and could estimate the duration 
of the atomic flights which is a time interval between two 
successive events when atom changes the sign of p. 

Just after SE at r = Tj, we have u — ► 0, z — > —1, 
p — > p + pj and H — ► Hj , and a change in the energy H 
during the time interval between the acts of SE is equal 
to the difference between the values of the pseudoenergy 
just before (Hj-\) and just after the j-th SE (Hj) 



Hj Hj. 



i= Hj - Hj 



-1 = UrP(Tj)pj + —Pj 
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+ u(Tj) COS x(Tj) + -jZ(Tj) + — (1 - Z )(Tj ~Tj-l), 

(5) 

where Hj is a value of the energy just after r = Tj 
and x(rj), u(tj-), z(rj), p(rj) are the values of the corre- 
sponding variables just before the moments r = Tj which 
are determined by coherent evolution between SE events. 
Changes in H at r = Tj are conditioned mainly by the cor- 
responding changes in u (Fig. 1). We stress that sudden 



changes in u occur at the moments of crossing the nodes 
of the standing wave and SE events. In Fig. lb the jump 
just after t = 1.2 /xs occurs when the atom crosses a node 
whereas the jump just before t — 1.4 (is is caused by a SE 
event. 

Let us estimate the average value of the energy jump's 
magnitude. Analytical estimate and simulation show that 
with sufficiently large values of the momentum. uj r \p\ > 
7/2, and at small detunings, we get (it cos x) ~ (z) ~ 
during the coherent evolution. The component u never 
goes far away from zero, and z performs frequency- 
modulated harmonic oscillations in the range — 1 < z < 1. 
The probability of SE is proportional to (z + l)/2. We es- 
timate the average value of the population inversion just 
before SE events at r = Tj to be equal to (z(rj)) ~ 0.5. In 
the expression (5) the first and fourth terms are estimated 
to be zero in average but the second, third, fifth and sixth 
ones are not. The total average change in Hj due to SE 
and the relaxation term is 



(Hj-H j . 1 ) = l ± + A ) 



(0) 



where we estimated the average value of z 2 to be 1/2 and 
the average value of the squared recoil momentum p 2 is 
1/3. The mean time between SE events is (Tj — Tj-i) = 
2/7. 

Thus, the evolution of the energy H is an asymmetric 
random walking. At positive values of the detuning A, 
the values of H increase in average, whereas at A < 
they may increase or decrease depending on the relations 
between uj t and A. We take the values of the detuning 
to be |A|/2 » oj r /9, and if A < then Hj decreases 
in average. The corresponding physical effects — light- 
induced acceleration and deceleration of atoms — are well 
known [1]. 

The diffusion coefficient in the energetic space can be 
estimated with the help of the two largest terms in Eq. 
(5) as 
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(7) 

where we suppose u to be a random-like process (Fig. 1) 
described by the Eq. (11) from Ref. [14], which is correct 
only under the condition \p\ <C l/uj r - This condition is 
fulfilled for all the atomic flights found in the numerical 
simulation. By the other hand, there is a stronger condi- 
tion, uj r \p\ > 7/2, that is needed to neglect correlations be- 
tween u 2 and cos 2 x. Thus, the expression (7) fails to give 
a result supporting the corresponding numerics for suffi- 
ciently slow atoms. In Fig. 2 we compare numerical and 
analytical results (obtained with Eq. (7) for the diffusion 
coefficient D with two values of the detuning A = —0.001 
and A = —0.0001. As expected, the correspondence is 
rather good for moderate values of the energy but not for 
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Fig. 2: The diffusion coefficient D in the energetic space vs 
the energy H: stars — A = -0.001, circles — A = -0.0001 
(numerical data) . Solid lines show the corresponding analytical 
curves obtained with the formula (7). 

small values (slow atoms) when the analytical formula (7) 
is not applicable. The drift velocity of particles in the 
energetic space at sufficiently large values of p is 

C " dr (rj-rj-x) 12 + 2 ' W 

This result is also correct only with sufficiently fast atoms. 
At the same condition as (8) we can estimate the friction 
force acting on the atom to be the following: 

dip) A7 , . 

F = -p. ~ — — . 9 
dr 2ui r p 

This kind of decreasing F with increasing p is a well-known 
fact [1]. 

Random jumps of the atomic energy just after SE give 
rise to a random walking of atoms in an optical lattice. In 
order to find distribution of the durations of atomic flights 
we consider the problem of the first passage time for the 
quantity Hj to return to its zero value (to be more correct, 
a return of H to \u\ must be considered, but with |A| <C 1 
the variable u cannot reach large values and always returns 
to zero value for the time ~ 2/7, see Fig. 1). In the very 
beginning of any flight we have Hj « 0, then it can reach 
a rather large value and after that it returns to zero. The 
duration of this process is a flight duration T. 

If the random jumps of Hj would be symmetric, the 
probability to find the flight duration to be equal to T 
would be proportional to T~ 15 , where the exponent does 
not depend on the diffusion coefficient (a classical result 
in theory of symmetric random walking [16]). 



What will happen if we take into account that the ran- 
dom walking of Hj is asymmetric? At A > 0, atoms begin 
to accelerate without any flights. At A < 0, the friction 
force tends to stop atoms, and instead of the power-law 
decay T~ 15 we get an exponential one at large T. At large 
times exceeding the mean SE time 2/7, one may treat the 
evolution of H as a diffusion process with a drift described 
by the Fokkcr-Planck equation in the energetic space 

dP d 2 P 

P{H , T) = - 2cm + D — m 

where the diffusion D and drift c coefficients are given by 
Eqs. (7) and (8), respectively. If the coefficients c and D 
would not change with changing the energy, the PDF for 
flight durations would be equal to [16] 

P fl oc e -o 2 T/D T -i.5 (n) 

This result agrees qualitatively with the results of numer- 
ical simulation shown in Fig. 3a for a few values of the 
detuning A. Really, all the PDFs in Fig. 3a have power- 
law fragments followed by exponential tails at large T in 
accordance with the formula (11). However, the length 
of these fragments depends strongly on the value of the 
detuning. At very small value A = — 1CP 5 (when coher- 
ent atomic dynamics is practically regular [14] and atom 
performs a random walk due to SE), Pa ~ T^ 15 , whereas 
at larger values of |A| the power-law fragments are much 
shorter. In fact, both c and D depend on the value of H, 
therefore, Eqs. (7) and (8) are not correct for small values 
of the momentum p, and more accurate formula for P^ is 
required. 

To illustrate the behavior of the atomic momentum at 
different values of the detuning we plot in Fig. 4a a typical 
chaotic walking of an atom with comparatively long flights 
at A = —0.001 and in Fig. 4b chaotic walking with short 
flights at A = -0.01. 

Qualitatively, the statistics in Fig. 3 are similar to those 
computed in Ref. [14] with a purely Hamiltonian coher- 
ent dynamics. Both of the PDFs contain power-law frag- 
ments with exponential decays at their tails, but the origin 
of these fragments and tails is different. Statistics of the 
purely Hamiltonian system [14] is governed by a determin- 
istic diffusion of the quantity u in a bounded space (be- 
tween its maximal and minimal values ±1), at a constant 
energy H . Solution for the first-passage-time problem for 
the quantity u = ±H gives a PDF with a power-law frag- 
ment with the slope —1.5 and an exponential tail. The 
quantity u jumps to zero value at the moments of SE (see 
Fig. 1), and its evolution cannot be treated as a diffusion 
process when we take into account SE. However, the en- 
ergy now is not a constant and can walk randomly (with 
a drift) within a broad range. The statistics of the system 
with SE is defined by a random walking of the energy H, 
not u. Thus, in both the systems the condition H as ±u 
defines the value of the energy which allows atoms to stop 
and turn back. Transformation of power-law fragments 
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Fig. 3: The PDFs Pa for the duration of atomic flights T ms in 
milliseconds with (a) small detunings (crosses A = — O.Of , stars 
A = -0.001, circles A = -0.0001, squares A = -0.00001) 
and (b) medium detunings (stars A = —0.09, a = —0.77; cir- 
cles A = -0.1, a = -0.27; squares A = -0.12, a = -0.05). 
Straight lines show slopes a of the power-law fragments of the 
PDFs in log-log scale. 



into exponential tails is explained in the purely Hamilto- 
nian system [14] by a limitation of the quantity u, whereas 
in the system with SE it is explained by a drift of the en- 
ergy H. Both of those factors prevent randomly walking 
quantities to go far from their critical values and decrease 
exponentially the probability of long atomic flights. 

Chaotic walking at moderate detunings, A| < 1. 

— The effect of velocity grouping, when there are one or 
a few values of the capture momentum p g to which current 
momenta of different atoms tend to, is known to occur at 
moderate negative values of the detuning [1]. This effect 
has been numerically demonstrated with Eqs. (3) in our 
recent paper [15]. Chaotic walking of atoms may occur 
in the regime of velocity groping if the values of p g are 
sufficiently small and atoms can change the direction of 
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Fig. 4: Dependencies of the current atomic velocity v on time: 
(a) A = —0.001, chaotic walking with comparatively long 
flights, (b) A = —0.01, chaotic walking with short flights, (c) 
A = —0.1, velocity grouping effect. 



motion due to fluctuations of momentum (caused by chaos 
in coherent evolution and/or SE, see Fig. 4c). 

Computed statistics of atomic flights at medium values 
of the detuning |A| < 1 (Fig. 3b) are similar to that at 
small detunings but the slope of the power-law decay may 
differ considerably. The PDFs P fl (T) shown in Fig. 3b 
demonstrate decrease of the slope with increasing the val- 
ues of | A with corresponding increase of the lengths of the 
power-law fragments. It should be emphasized that the 
length of the power-law fragments (and the mean flight 
length) increases significantly with a rather small increase 
in |A|. 

Varying the value of the detuning A, one can manipu- 
late atomic transport in an optical lattice and its statis- 
tical properties. In Fig. 5 we plot the dependencies of 
the mean duration of atomic flights (T) (in microseconds) 
and the slope of the power-law fragments of Pr ~ T a on 
the detuning in the range of medium values. The depen- 
dencies correlate well with each other. Figure 3 b demon- 
strates clearly that the length of the power-law fragments 
increases with increasing the absolute values of the detun- 
ing, whereas the absolute value of the slope a decreases. 
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Fig. 5: Dependencies of the logarithms of the mean duration 
of atomic flights (T MS ) (in microseconds) (solid line) and of 
the slope a of the PDF power-law fragments (squares) on the 
medium detuning A. 



The mean duration of flights (T) increases correspondingly 
with increasing |A|. 

The control is nonlinear in the sense that when slightly 
decreasing A, say, from —0.08 to —0.12, the mean time of 
flights increases in a few orders of magnitude (see Fig. 5). 
This effect is explained by increase in the value of the 
capture momentum p g and decreasing fluctuations of the 
current atomic momentum with increase in the absolute 
value of |A|. For example, at A = —0.1, one gets p g ~ 500 
(~ 1.5 m/s) with the momentum fluctuations of the same 
order. Thus, the atom can change its direction of motion 
(see Fig. 4 c, with A = —0.1 the mean time of flights 
is T ps ~ 10/xs) but not so frequently as in the case of 
smaller detunings (see Fig. 4b). The value of p g increases 
with increasing |A| and, say, at p g ~ 1000 the momentum 
fluctuations are of the order of 300, and the atom cannot 
change its direction of motion. Reduction in the momen- 
tum fluctuations is caused by increase in the friction force 
F. With the values of p smaller than p g , the force F is 
so large that changes in the energy H with time is not a 
process of random walking (as it is in the case with smaller 
detunings) but rather a directed drift in the momentum 
space to the value of the capture momentum. Thus, it is 
practically impossible for atoms to decrease their current 
values of p to zero value, and the process of chaotic walk- 
ing eventually stops. The slope a of the PDF power-law 
fragments can go to zero due to existence of an exponen- 
tial decay at the very tail of Pa- With a purely power-law 
decay the minimal slope would be a = — 1. 

In conclusion, we have shown that near the resonance 
atomic transport in an optical lattice is a complicated pro- 
cess of chaotic walking caused by an interplay between 



coherent but deterministically chaotic internal atomic dy- 
namics and spontaneous emission random events. It is 
possible to manipulate this process and its statistical char- 
acteristics by varying the single control parameter, the 
atom-laser detuning A. 
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